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^ ! Abstract 

.^ ' The effects of removing large scales external to the inertial range on the prop- 

^^ . erties of scales within the inertial range are studied in a high-Reynolds-number 

p^ I turbulent flow. Structure functions of both even and odd orders are strongly af- 

Q>^ ' fected across the entire inertial range, but odd-order moments are affected to a 

0^ ■ greater degree. In particular, the skewness of velocity increments shows a signifi- 

^ , cant reduction whereas the flatness changes comparatively little. The reduction in 

t^' skewness is counterbalanced essentially by the interaction between the small-scale 

I ' energy and the large-scale rate of strain. The implications of these results for the 

O , conventional cascade picture are examined briefly. 

^ : 
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d 1 Introduction 



The Karman-Howarth equation^ is an exact dynamical equation for the second- 
order correlation function in isotropic turbulence. The equation is easily expressed 
in terms of the longitudinal structure function S2{r,t), where 5*2 = (An^) = {[u{x + 
r) — n(x)]^), u{x) is the turbulent velocity component in the direction x, and r is 
the separation distance in that same direction. Assuming that the time variation of 
S2 is slow, Kolmogorov^ reduced the Karman-Howarth equation to an equation for 
the third-order structure function. This equation, given by 

Ssir) ^ {[u{x + r) - u{x)f) = -^{e)r + 61^^^, (1) 



has been rederived under less restrictive circumstances^''^, and works well in shear 
flows as long as r is much smaller than a characteristic large scale, L. Here {e) is 
the average value of the energy dissipation rate. Figure 1 is a typical comparison 
between the measured third-order structure function and the right hand side of 
Eq. (1). In particular, there is a range of scales over which the viscous term is 
negligibly small, and we have 

S,{r) = --^{e)r. (2) 

Indeed, it is conventional to define the inertial range as the range of scales for which 
this equation is valid. 

The usual interpretation of Eq. (2) is that, on the average, it prescribes a uni- 
directional down-scale energy flux.^'^ It is consequently thought that the statistical 
physics appropriate to the inertial range has a strongly non-equilibrium character.^ 
On the other hand, consider the following argument. If the inertial range dynamics 
is not influenced by viscosity (as apparent in the previous paragraph), it is rea- 
sonable to suppose that the dynamics there is governed essentially by the Euler 
equations. The Euler equations possess time-reversal symmetry, t — > — t, u -^ — u 
(t = time, u = vector velocity), which could therefore be interpreted to mean that 
the energy flow in the inertial range is essentially as much up-scale as down-scale. 
The proper meaning of Eq. (2) would then be that the energy flux it prescribes is a 
small difference between equally large energy flows in the up-scale and down-scale 
directions. The finiteness of 5*3 may then mean that there exists no more than a 
"small leakage" between two opposing and equally strong fluxes.'^ If this is true, the 
situation calls for a simpler type of statistical physics. 

Prompted by such qualitative considerations, we studied, some nine years ago,^ 
the source of energy flux in the Kolmogorov equation. A result of particular interest 
is the experimental finding that ^3 became small when large scales in the fluctuating 
velocity were filtered out. This result was thought to imply that the Kolmogorov 
result, Eq. (2), was more a reflection on the source of the energy fiux (large scales) 
than a statement about unidirectional cascades. We present here an up-date of the 
experimental part of that (unpublished) study. We use velocity signals acquired in a 
very high-Reynolds-number atmospheric surface layer, described below, to analyse 
the effect of removing large scales on the statistics of An,., in particular the effects 
on its third moment. The principal conclusion is that, while the removal of large 



scales affects both even-order and odd-order moments of velocity increments in 
the inertial range, the odd moments are much more vulnerable. The results are 
interpreted briefly in the light of the issues raised above, though it is fair to warn 
the reader that we shall not be able to resolve the issue satisfactorily. Section II 
is a brief commentary of the experimental data, while Sec. Ill presents the basic 
experimental results. The chief conclusions are discussed in Sec. IV. 

2 The data and the filter characteristics 

The filtering experiments to be described below have been made for a number 
of flows such as a pipe flow (diameter Reynolds number ~ 230,000), boundary 
layer turbulence (boundary layer thickness Reynolds number ~ 20,000), atmospheric 
turbulence at a height of 6 m above the ground (Taylor microscale Reynolds number 
R\ ~ 2,000), and atmospheric turbulence at a height of 35 m above the ground 
{Rx ~ 10,000).^ The results are entirely consistent with each other, but are more 
revealing when the Reynolds number is high and the scale separation large, so we 
present data only for the last case. Tests were made on several sets of data, but the 
results are presented for brevity for only one. 

The atmospheric data were acquired by a hot-wire probe mounted at a height 
of 35 m above the ground on the meteorological tower at the Brookhaven National 
Laboratory. The wind speed and direction were independently monitored by a 
vane anemometer mounted close to the tower. During data acquisition, the mean 
wind was essentially steady in both magnitude and direction. This was ensured, 
a posteriori, by using only those sets of data for which the averages computed 
over different segments of the measurement interval were steady. The conditions of 
measurement thus approached 'controlled' circumstances. 

The hotwire, about 0.7 mm in length and 5 fim in diameter, was operated in a 
constant temperature mode. Its frequency response was compensated to be flat up 
to 20 kHz. The voltage from the anemometer was low-pass filtered and digitized. 
The low-pass cut-off was half the sampling frequency, fs, which was set at 5,000 Hz. 
This frequency was high enough to include essentially all the finest-scale turbulent 
fluctuations in the flow. The voltage was constantly monitored on an oscilloscope to 
ensure that it did not exceed the digitizer limits. The voltage from the anemometer 
was converted to wind velocity through the standard calibration procedure, which 



included calibrating the hotwire just prior to mounting and checking it immediately 
after dismounting. Taylor's frozen flow assumption was employed to convert the 
time variable into the spatial distance. 

The following characteristics pertain to the present data: The mean wind speed 
U ~ 7.6 ms^-*^, root-mean-square fluctuation velocity u' ~ 1.36 ms^^, the mean 
energy dissipation rate (e) ^ 3 x 10^^ m^s~^, the Kolmogorov scale rj ~ 0.57 
mm, the Taylor microscale in the streamwise direction A ~ 11.4 cm. The Taylor 
microscale Reynolds number u'X/i^ ~ 10, 340. We use the height from the ground, 
H, as a measure of the large scale. If we assume that the classical boundary layer 
arguments are valid, a characteristic scale for the large-scale transport is of the order 
0.4 X 35 m, where 0.4 is the 'accepted' numerical value of the Karman constant. 

To obtain the filtered signal we used a Butterworth filter with excellent cut-off 
characteristics and no phase shift. The filtering operation was done by convolution 
in real space. The cut-off scale rj and the cut-off frequency fo are related by 
Vf = U/fo- After filtering in the forward direction, the filtered sequence is reversed 
and run back through the filter. The resulting sequence has no phase distortion and 
double the filter order. We have used filters of order 4 and 10 (i.e., of order 2 and 5, 
respectively, but with two passes), mostly the latter. The filter characteristics are 
demonstrated in Figs. 2 to 4. (A different filtering scheme was employed in Ref. [8] 
and the nature of the results was identical.) Figure 2 shows the energy spectral 
density of streamwise velocity fluctuations. The Kolmogorov scaling of the spectral 
density for the unfiltered signal extends roughly over three decades of scales. The 
spectrum for the filtered signal, with high-pass setting at 2 Hz is also shown for the 
tenth-order filter. The attenuation in Fourier space exceeds 40 dB per octave. The 
ratio of the power spectral density of the filtered signal to that of the unfiltered 
signal is shown in Fig. 3 for fourth-order and tenth-order filters. For the latter, the 
filter effects vanish by within a factor of 2 of the cut-off setting. 

The filter produces minimal phase distortion, as demonstrated below. Let u be 
the unfiltered velocity, and u-^ and n^ be the high-pass (small-scale) and low-pass 
(coarse-scale) parts of u, both obtained by independent filtering operations with the 
filter set at fo Hz. Define uf=u — u^. If the filter does not produce any phase 
distortion, we should find u^ = uf. Figure 4 shows this is indeed so. 



3 Experimental results 

3.1 Probability density functions 

It is known^'^"-'^^ that the probabihty density function (pdf) of the velocity incre- 
ment, Aur, continually changes its character as the separation distance r is varied: 
it is roughly square-root exponential for r in the dissipative range and Gaussian 
when r is comparable to the large scale. The pdfs are shown in Fig. 5 (a). Figures 5 
(b), (c) and (d) show similar data for high-pass filtered signals, each for a different 
filter setting. The spatial cut-off scale ry is given in the caption. The removal of the 
large scales changes the pdfs in two respects: (1) the shape is affected significantly 
for separation distances r smaller than the smallest scale removed by filtering, the 
latter being of the order rj; (2) the asymptotic form of the pdf, for r > rf, does 
not differ much from that at r = rj, and is far from Gaussian. These effects can 
be quantified by fitting a stretched exponential, pAur = p{0)exp[—c\Aur\"^^^'], to 
the pdfs of filtered signals, and comparing the stretching exponent m(r) for var- 
ious filter settings. This is done in Fig. 6. The exponent m is slightly different 
for positive and negative sides of the distribution (shown by full and dashed lines, 
respectively), but these differences do not mask the major effect: the exponent is 
noticeably smaller for filtered signals than for the unfiltered signal, and the effect 
occurs at scales significantly smaller than the smallest filtered scale. 

3.2 The third-order structure function 

Figure 7 is a summary of the effects of the filtering operation on the third-order 
structure function. The ordinate is normalized by r(e), and the separation distance 
is normalized by the Kolmogorov scale rj for the unfiltered data. First, a remark is 
in order on the function K = (Au^)/r(e) for the unfiltered data: the inertial range, 
defined as the region where K is flat and equal to 4/5, ranges roughly between lOOr/ 
and 2500rj. (Note that the spectral density shown in Fig. 1 appears to possess a 
larger scaling range, but it is more prudent to base this estimate on the basis of 
Eq. (2).) 

The high-pass filtered data are obtained by removing all scales above 30 m 
(frquencies below about 0.25 Hz), 3.8 m (~ 2 Hz) and 1.5 m (--^ 5 Hz). The arrows 
mark the values of rj/77. For the 3.8 m case, for example, (Anj!) is smaller by an 
order of magnitude for r/rf = 0.5, by a factor 2 for r/rf = 10, and by about 10% 



for r/rj = 100. Inferring roughly from Fig. 2 that the artifacts of the filtering 
operation do not extend to scales smaller than about rj/2, it is clear that the 
removal of large scales beyond the upper edge of the inertial range (~ 1.7 m) has a 
strong effect throughout the inertial range. However, the dissipation scales are not 
affected significantly, so (e) and ij remain essentially the same as for the unfiltered 
data; we have therefore simply used the unfiltered values for normalization. The 
dissipation scales just begin to feel the effect of filtering for rj = 1.5 m, but this 
cut-off begins to encroach directly on the inertial range itself. For most of the results 
below, we shall use rj = 3.8 m (--^ /o = 2 Hz), unless stated otherwise. 

A naive interpretation of Fig. 2 is that the large scales contribute significantly to 
(Anj!) in the inertial range. A closer look demands a more complex interpretation. 
By splitting, as before, the velocity signal u into the fine-scale (high-frequency) 
and coarse-scale (low-frequency) velocity components u-^ and n^, separated by the 
cut-off scale rj or the cut-off frequency fo, and noting that u = u-^ + u"^ and 
Aur = Au^ + Au^, the third-order structure function can be written as 

{A4) = {iAu>f) + 3{iAu>fAu<) + 3{iAu<fAu>) + {iAu<f). (3) 

In the above equation, it is not hard to argue that ((An^)^) and 3((An^)^Au^) 
should both be small for very small r, and that the term ((An^)^) should approach 
(Anj!) itself as r ^ 0. For intermediate values of r, it is not easy to estimate 
the correct orders of magnitude of these terms. Experimental data (see Fig. 8) 
show that the largest term is ((Au^)^) in the lower part of the inertial range and 
3((An^)^Au3 in the upper part, until they are both overtaken, as r — > rj, by the 
exclusively large-scale part ((Au^)^). 

The dominance of the term 3((An^)^Au3 suggests that a significant part of 
the third-order structure function comes from the straining of the small-scale kinetic 
energy, (An^)^, by the large-scale feature Au^. This suggests the existence of a 
coupling between scales that are significantly separated. The effect is illustrated in a 
different manner in Fig. 9 where the second order moment of Aur, conditioned on the 
velocity gradient of the cut-off scale, du^/drf, is plotted against the latter quantity; 
the conditioning variable is the gradient of the largest scale rejected by the filtering 
operation. The variance may be expected to be independent of the conditioning 
variable if the energy transfer across scales were purely a local phenomenon — in 
contrast to the measured effect. 
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3.3 Odd versus even moments 

The second and fourth order structure functions for the unfiltered and filtered data 
are shown in Fig. 10. Filtering has some effects on even-order moments as well; the 
larger the moment order, the larger this influence. However, these effects are much 
weaker than that for the third-order moment (and for other odd-order moments as 
well, though not shown here). Two specific statements can be made. First, for the 
second-order structure function, the filter is felt only up to scales that are at best 
an order of magnitude smaller than r j , but the effect on the third-order extends to 
scales that are smaller by one further order of magnitude. Second, the effects on 
the third-order are larger than those on the next-higher even moment, namely the 
fourth. 

A better appreciation for differences between odd and even moments can be 
gained by examining normalized moments. Figure 11 shows the skewness of the 
velocity increment for the three filter settings used above, as well as for the un- 
filtered data. A comment first on the unfiltered signal itself: the fact that the 
skewness of the unfiltered signal is not constant in the inertial range shows that the 
r-dependence of (Au^) has a power-law exponent that is higher than 2/3. While a 
rough estimate of this "anomalous correction" for second-order structure function 
can be obtained from this graph, it has been argued elsewhere^ that a proper esti- 
mate needs greater care. This is not our main point, so we shall make no further 
reference to it here. Note that the effect of removing scales above 30 m (below 0.25 
Hz) already penetrates the inertial range; removing scales above 3.8 m (below 2 
Hz) and 1.5 m (5 Hz) significantly reduces the skewness of scales everywhere in the 
inertial range. For instance, the skewness for the 2 Hz case is reduced by a factor of 
2 for scales one order of magnitude smaller than rj. On the other hand, the fiatness 
of Aur shows hardly any effect for scales smaller than rf/2 (Fig. 12). 

4 Discussion and conclusion 

A principal conclusion of this work is that the removal of large scales significantly 
affects the statistics well into the inertial range. In terms of the normalized moments 
of velocity increments, the skewness — as are other odd moments — is affected much 
more strongly than the flatness — or other even moments. Indeed, if the large scales 
can be thought of as related, loosely, to the "size of the system" , it appears that the 



system-size effects are felt more strongly on odd moments than on even moments 
of velocity increments. 

It is worth asking if this effect is special to shear flows, because all the flows we 
have explored are shear flows. But we have some experience in assessing the shear 
effects^ on scaling, and this experience suggests that the present conclusion is not 
restricted. It would, however, be worth repeating these calculations for isotropic 
turbulence. We have not done so because all the isotropic data available to us — 
experimentally as well as by simulations — possess limited scale separation, which 
makes i t very difficult to extract meaningful results. 

It must be emphasized that the present results do not suggest that the skewness 
of the velocity increment resides in the large r. Rather, what they suggest is that this 
asymmetry of the pdf is related to the straining of the energy in the inertial range 
by large scales. The asymmetry is the subject of another study in its own right, ^^ 
but it should be pointed out that it is associated with all magnitudes of An,.. To see 
this, consider the difference between the two sides of the pdf oi z = Aur/{r{e))^'^. 
The quantity z^ x [p{z) — p{—z)], when integrated over < z < oo, yields the 
function, K = {Auf)/r{£). This integrand z^ x [p{z) —p{—z)] is plotted against z in 
Fig. 13. The figure shows, for these specific conditions, that most of the contribution 
to K comes largely from z lying between about 2 and 10. The small values of z 
lying between and 2 make the opposite contribution to K. The filtering operation 
essentially affects the integrand everywhere. 

What do the present results imply for the energy cascade and the nature of non- 
equilibrium (or otherwise) prevailing in the inertial range? They show that there are 
significant nonlocal effects, and that the energy cascade, even as an average concept, 
is not driven entirely by local effects. While cascade models do have pedagogi 
cal value, they seem to miss this element. However, interpretations on the non- 
equilibrium nature of inertial-range dynamics are difficult because the information at 
hand is very limited. While it is unclear that the time reversal symmetry of the Euler 
equations guarantees the near-equality of fluxes in both directions (for instance, 
Euler equations can in principle produce energy flux by an inviscid mechanism^^), 
it is equally unclear that Eq. (2) in itself demands a strong non-equilibrium situation 
in the inertial range. 

As a contrast, it is worth remembering that the subtle manner in which irre- 
versibility appears even in the simpler case of the kinetic theory of gases took a long 



time to unfold. One should therefore not feel too pessimistic about the inconclu- 
siveness of the present situation. 
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Figure 1: The (negative value of tlie) third-order structure function, 6*3, plotted against the 
separation distance, r, normalized by the Kolmogorov scale r/. The measured data are compared 
with Eq. (1), full line. The agreement is quite good up to r = iJ/20. The largest disagreement, 
seen at the very small scales, is due to the difficulty of resolving those scales well. The dashed 
line, which differs from Eq. (1) only below rJT] of about 50, is Eq. (2). The data are for the 
atmospheric surface layer at the height iJ = 35 m above the ground, and all averages are over 
time; details of measurements are given in Sec. II. Here, H is treated as the representative large 
scale, and marked by an arrow on the abscissae. 
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Figure 2: The full line is the power spectral density (PSD) of the longitudinal (streaniwise) 
velocity fluctuation. The dashed line is the effect of a high-pass filter that is set to /o = 2 Hz. 
The filter is sharp and does not produce oscillations. 
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Figure 3: The ratio of the filtered to the unfiltered power spectral densities for the fourth and 
tenth order filters. We have experimented with both filters but used the tenth-order filter for 
the data presented in the rest of the text. For this filter, the filter effects essentially disappear 
for frequencies smaller than about /o/2. 
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an independent filtering operation, also performed with fo set to 2 Hz (heavy dashed line). The 
agreement between the two is excellent, as should be expected if the filtering operation produced 
no significant phase distortion. 
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Figure 6: The stretching exponent, m, versus the separation distance for the unfiltered and fil- 
tered data sets. The values of m are obtained by fitting empirically pam^ = p(0)exp[— cjAurl™^'"-'] 
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Figure 7: The Kolmogorov function (Au^)/r(e) plotted for unfiltered ( — ) and high-pass filtered 
traces (x, /„ = 0.25 Hz, rj/H Ki 0.88; o, /„ = 2 Hz, ry/if Ri 0.11; □, /o = 5 Hz, r//F « 0.044). 
The normalizing scales r] and (e) are for the unfiltered signal, but the filtering does not change 
them to any significant extent. The horizontal line extending between r/ij of 100 and 2500 is 
the estimated inertial range. In this figure and others to follow, the cut-off scale rj is shown 
by an arrow on the abscissae. The removal of large scales affects scales that are typically two 
decades smaller than rf. 
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Figure 8: The result of decomposing the third-order structure function (An^) in terms of Au^ 
and Au>; /^ = 2 Hz, r//i7 « 0.11. o, {Auf)/r{e); □, ((Au>)3)/r(e); +, 3{{Au>fAu<)/r{e); 
A, 3((Au<)^Au^)/r(e); o, ((Aii^)'^)/r(e). In the upper part of the inertial range, the most 
dominant term is 3((Au^)^A«^ ra. 
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Figure 9: The conditional expectation of the unfiltered variance (Au^.)^ conditioned on the rate 
of strain of the largest scale retained in the filtering operation, dAu^/dx; fo = 2 Hz, rf/ij ~ 
7,000. The various symbols correspond to different scale separations r, as noted in the inset. 



19 




Figure 10: Second and fourth order structure functions for tlie unfiltered (lines) and filtered 
data (symbols). — and x are for the second-order; — — and o are for the fourth-order; /o = 2 
Hz, Tf/ri K, 7,000. 
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Figure 11: The skewness of velocity increments for the unfiltered ( — ) and filtered signals 
(x,/o = 0.25 Hz, Vf/H « 0.88; o,/^ = 2 Hz, rf/H w 0.11; □, /„ = 5 Hz, rj/ilf « 0.044). The 
filter settings are shown by arrows. 
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Figure 12: The flatness factor of velocity increments for the unflltered ( — ) and filtered signals 
(x,/o = 0.25 Hz, Vf/H « 0.88; o,/^ = 2 Hz, rj/H « 0.11; □, /o = 5 Hz, rj/H « 0.044). The 
filter settings are shown by arrows. 
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Z = Au/(r<e))' 

Figure 13: The integrand of K = {Aur)^/r{e), plotted against Aur/{r{e))^'^. The crosses 
are for the unfiltered signal, r/r] = 1000, and — represents a smooth line through them. The 
corresponding quantity for the filtered data, with fo = 2 Hz, rf/r] ~ 7000, are shown by 
diamonds. The effect of the filter is drastic at all scales. The relatively large scatter in the data 
is the result of taking the difference between two nearly equal quantities. 
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